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SCHEMES FOR COMPUTING PERFORMANCE PARAMETERS OF 
FIBER OPTIC GYROSCOPES 

CROSS REFERENCE TO RELATED APPLICATION 

This application claims the priority of U.S. provisional Patent Application 

5 Serial No. 60/442,634 (by Humphrey, filed January 24, 2003, and entitled 
"SCHEMES FOR COMPUTING PERFORMANCE PARAMETERS OF FIBER 
OPTIC GYROSCOPES"). 
Background 

(1) Field 

10 The present disclosure relates to schemes for computing performance 

parameters of fiber optic gyroscopes (FOGs) using closed-loop transfer functions. 

(2) Description of Related Art 

A FOG is a device that can detect rotation in a variety of applications, 
including navigation and stabilization schemes. Generally, a FOG can include an 

15 optical subsystem and an electrical subsystem. The optical and electrical 
subsystems can provide inputs to each other. 

A FOG can be characterized by a variety of performance parameters, 
including an operating frequency and a bandwidth. Generally, schemes for 
computing FOG performance parameters separately model FOG optical and 

20 electrical subsystems with two open-loop systems. Since FOGs can operate with 
their optical and electrical subsystems in a closed-loop configuration, however, 
meaningful conclusions cannot be reliably provided by two open-loop systems. 
Summary 

Schemes for computing performance parameters of FOGs using closed- 
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loop transfer functions are described herein. 

A method for computing a performance parameter of a FOG is described 
herein. In one embodiment, the method may include providing a closed-loop 
transfer function based on optical components and electrical components of the 
5 FOG; based on the transfer function, determining a relationship between the 
performance parameter and at least one physical parameter associated with at least 
one component of the FOG; and, based on the relationship, computing the 
performance parameter. 

In one aspect, providing may include providing a feedforward component 
10 representing at least one FOG optical component and at least one FOG electronic 
component; and, providing a feedback component representing at least one FOG 
optical component and at least one FOG electronics component. 

In one aspect, providing a feedforward component may include 
representing, in the feedforward component, at least one noise component. 
15 In one aspect, providing a feedforward component may include 

representing, in the feedforward component, at least one disturbance, wherein the 
at least one disturbance is based on at least one of: an optical power noise, a shot 
noise, a preamplifier current noise, a preamplifier thermal noise, a preamplifier 
voltage noise, and an analog-to-digital converter (ADC) quantization noise. 
20 In one aspect, providing a feedforward component may include 

representing, in the feedforward component, at least one of: a phase modulator, a 
photodetector and an associated preamplifier, a filter, an ADC, and a sampler. 

In one aspect, representing the phase modulator may include representing 
the phase modulator based on an optical power of a light beam propagating 
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through a fiber-optic coil and an operating phase bias. 

In one aspect, representing the phase modulator may include representing 
the phase modulator based on a product of the optical power and a sinusoidal 
function of the operating phase bias. 
5 In one aspect, representing the photodetector and the associated 

preamplifier may include representing the photodetector and the associated 
preamplifier based on a photodetector scale factor, a preamplifier impedance, and 
a preamplifier gain. 

In one aspect, representing the photodetector and the associated 
10 preamplifier may include representing the photodetector and the associated 

preamplifier based on a product of the photodetector scale factor, the preamplifier 
impedance, and the preamplifier gain. 

In one aspect, representing the filter may include representing the filter as a 
gain in voltage after the photodetector and associated preamplifier and before the 
15 ADC 

In one aspect, representing the ADC may include representing the ADC as 
a gain based on the number of bits in the ADC. 

In one aspect, providing a feedback component may include representing, 
in the feedback component, at least one of: sampler, a truncator, a digital-to- 
20 analog converter (DAC), a phase modulator, and a fiber-optic coil. 

In one aspect, representing the truncator may include representing the 
truncator as a digital truncation gain. 

In one aspect, representing the DAC may include representing the DAC as 
a gain based on the number of bits in the DAC. 
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In one aspect, representing the phase modulator may include representing 
the phase modulator as a scale factor. 

In one aspect, representing the fiber-optic coil may include representing the 
fiber-optic coil as a time delay. 
5 In one aspect, representing the fiber-optic coil may include representing the 

fiber-optic coil based on a transit time for a light beam to propagate through the 
fiber-optic coil. 

In one aspect, determining a relationship may include, based on the transfer 
function, determining a relationship between the performance parameter and at 
10 least one physical parameter associated with at least one component of the FOG, 
wherein the at least one physical parameter includes at least one of: 

an optical power of a light beam propagating through a fiber-optic coil, an 
operating phase bias, a photodetector scale factor, a preamplifier impedance, a 
preamplifier gain, a filter gain, an ADC gain, a digital truncation gain, a DAC 
15 gain, a transit time for a light beam to propagate through the fiber-optic coil, and a 
phase modulator scale factor. 

In one aspect, computing may include providing an input based on a rate of 
rotation of a fiber-optic coil and a scale factor, the scale factor including a 
wavelength of a light beam propagating through the coil, a coil length, and a coil 
20 diameter. 

In one aspect, computing may include computing a performance parameter 
including at least one of a bandwidth, a coefficient of random walk, an operating 
frequency, and a power spectral density of noise. 

In one embodiment, the method may further include providing a value of a 
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performance parameter and determining at least one value associated with the at 
least one physical parameter for which the computed performance parameter will 
have the value. 

In one aspect, determining the at least one value may include providing at 
5 least one initial value associated with the at least one physical parameter; based on 
the relationship and the at least one initial value, computing the performance 
parameter; and, based on a difference between the computed performance 
parameter and the value, iteratively adjusting at least one value associated with the 
at least one physical parameter and iteratively computing the performance 
10 parameter. 

In one embodiment, the method may further include providing a first value 
of a first performance parameter; providing a second value of a second 
performance parameter; and, determining at least one value associated with the at 
least one physical parameter for which the computed first performance parameter 

15 will approach the first value and the computed second performance parameter will 
approach the second value. 

In one aspect, determining at least one value may include providing at least 
one initial value associated with the at least one physical parameter; based on the 
corresponding relationship and the at least one initial value, computing the first 

20 performance parameter and the second performance parameter; and, based on a 
difference between at least one of the first value and the computed first 
performance parameter and the second value and the computed second 
performance parameter, iteratively adjusting at least one value associated with the 
at least one physical parameter and iteratively computing the first performance 
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parameter and the second performance parameter. 

A processor program for computing a performance parameter of a fiber- 
optic gyroscope (FOG) is described herein. In one embodiment, the processor 
program may be stored on a processor-readable medium and may include 
5 instructions to cause a processor to receive a closed-loop transfer function based 
on optical components and electrical components of the FOG; based on the 
transfer function, determine a relationship between the performance parameter and 
at least one physical parameter associated with at least one component of the 
FOG; and, based on the relationship, computing the performance parameter. 
10 In one aspect, the instructions to compute may include instructions to 

compute a performance parameter including at least one of a bandwidth, a 
coefficient of random walk, an operating frequency, and a power spectral density 
of noise. 

In one embodiment, the processor program may also include instructions to 
15 receive a value of a performance parameter, and determine at least one value 

associated with the at least one physical parameter for which the computed 

performance parameter will have the value. 

In one aspect, the instructions to determine may include instructions to 

receive at least one initial value associated with the at least one physical 
20 parameter; based on the relationship and the at least one initial value, compute the 

performance parameter; and, based on a difference between the computed 

performance parameter and the value, iteratively adjust at least one value 

associated with the at least one physical parameter and iteratively compute the 

performance parameter. 
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In one embodiment, the processor program may also include instructions to 
receive a first value of a first performance parameter; receive a second value of a 
second performance parameter; and, determine at least one value associated with 
the at least one physical parameter for which the computed first performance 
5 parameter will approach the first value and the computed second performance 
parameter will approach the second value. 

In one aspect, the instructions to determine may include instructions to 
receive at least one initial value associated with the at least one physical 
parameter; based on the corresponding relationship and the at least one initial 

10 value, compute the first performance parameter and the second performance 
parameter; and, based on a difference between at least one of the first value and 
the computed first performance parameter, and the second value and the computed 
second performance parameter, iteratively adjust at least one value associated with 
the at least one physical parameter and iteratively compute the first performance 

15 parameter and the second performance parameter. 

These and other features of the schemes for computing performance 
parameters of FOGs described herein may be more fully understood by referring 
to the following detailed description and accompanying drawings. 
Brief Description of the Drawings 

20 Fig. 1 is a block diagram of an exemplary closed-loop transfer function for 

a FOG. 

Fig. 2 is a block diagram of an exemplary feedforward component of the 
closed-loop transfer function shown in Fig. 2 

Fig. 3 schematically illustrates a prior-art FOG. 


8 NGC- 162/ 000388-280 

Detailed Description 

Certain exemplary embodiments will now be described to provide an 
overall understanding of the schemes for computing performance parameters of 
FOGs described herein. One or more examples of the exemplary embodiments are 
5 shown in the drawings. 

Those of ordinary skill in the art will understand that the schemes for computing 
performance parameters of FOGS described herein can be adapted and modified to 
provide devices, methods, schemes, and systems for other applications, and that 
other additions and modifications can be made to the schemes described herein 

10 without departing from the scope of the present disclosure. For example, 

components, features, modules, and/or aspects of the exemplary embodiments can 
be combined, separated, interchanged, and/or rearranged to generate other 
embodiments. Such modifications and variations are intended to be included 
within the scope of the present disclosure. 

15 Generally, the exemplary schemes described herein include a closed-loop 

representation of FOG optical subsystem components and FOG electrical 
subsystem components to compute performance parameters for FOGs. In one 
embodiment* a closed-loop transfer function can be used to determine a 
relationship between a FOG performance parameter and physical parameter(s) 

20 associated with FOG component(s). The relationship can be used to determine 
value(s) of the physical parameter(s) for which the performance parameter will 
approach a performance parameter value. 

Fig. 3 schematically illustrates a prior-art FOG. FOGs are well known and 
may be understood by referring to the disclosures of U.S. Patent Nos. 4,705,399 to 
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Graindorge et al. and 5,337,142 to Lefevre et al., the contents of which patents are 
expressly incorporated by reference herein. 

As shown in Fig. 3, FOG 10 may include an optical subsystem 12 and an 
electrical subsystem 14. Optical subsystem 12 may include a light source 22, a 
5 beam splitter 24, a phase modulator 26, and an optical waveguide 28. Electrical 
subsystem 14 may include a signal digitizer 30 and a demodulator 32. Optical 
subsystem 12 can provide a signal 16 to electrical subsystem 14, and electrical 
subsystem 14 can provide a feedback signal 18 to optical subsystem 12. Electrical 
subsystem 14 can also provide a signal 20 to an application. FOG components 22, 
10 24, 26, 28, 30, and 32 may be connected by optical and/or electrical connection(s) 
and may communicate with component(s) other than those illustrated. 

Operation of FOG 10 may be briefly understood in the following manner. 
Light source 22 can provide a light signal 15 to beam splitter 24, and beam splitter 
24 can split the light signal into two light signals that travel in opposite directions 
15 34, 36 along an optical path defined by optical waveguide 28. Beam splitter 24 can 
receive the two light signals exiting from optical waveguide 28, combine the two 
light signals, and provide the combined light signal 16 to signal digitizer 30. Based 
on the combined light signal 16, signal digitizer 30 can produce an output signal 
proportional to a phase difference between the two light signals exiting the optical 
20 waveguide 28. According to the well known Sagnac effect, this phase difference 
can be used to measure a rate of rotation of the optical waveguide 28. 

A variety of schemes for adjusting the operating point of a FOG 10 are 
available. Generally, these schemes superimpose artificial phase differences on 
the two light signals 34, 36 counterpropagating in the optical waveguide 28. In 
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these schemes, the output from the signal digitizer 30 can be provided to the 
demodulator 32, and the demodulator 32 can provide a feedback signal 18 to phase 
modulator 26 to modulate the relative phases of the counterpropagating light 
beams. 

5 Fig. 1 is a block diagram of an exemplary closed-loop transfer function for 

FOG 

10. As shown in Fig. 1, the transfer function 100 may include an input 110, a 
summing point 120, a feedforward component 130, a feedback component 140, 
and a branch point 150. Input 110 and feedback component 140 may be provided 

10 to positive and negative terminals 122, 124 of summing point 120, respectively. 
As described herein, transfer function 100 may be used to compute an operating 
frequency and bandwidth of FOG 10. Appendices 1-5 include features of transfer 
function 100 described herein. 

Generally, input 110 may be based on a rate of rotation of an optical 

15 waveguide 28 and a scale factor. Input 110 may be based on a product or the rate 
of rotation and the scale factor. In one embodiment, the scale factor may include a 
wavelength of light propagating through the optical waveguide 28, an optical path 
length of the optical waveguide 28, and a diameter of the optical waveguide 28. 
The scale factor may be associated with the well known Sagnac scale factor. For 

20 example, in one embodiment of FOG 10, an optical waveguide 28 may include a 
coil of optical fiber wound on a spool-type structure, such as a bobbin, and a light 
source 22 that can be, for example, a superluminescent diode (SLD). In such an 
embodiment, the input 110 may be represented as the product 

(Eq. 1) GK.=<^ 
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where Q is the rate of rotation of the coil, K s is the well known Sagnac scale 
factor, L is the length of the coil, D is the diameter of the coil, X is the wavelength 
of light emitted by the SLD, and c is the speed of light in vacuo. 
5 Feedforward component 130 may include representations of at least one FOG 
optical component and at least one FOG electrical component. As shown in Fig. 1, 
feedforward component 130 may include a representation 132 of a phase 
modulator 26. In one embodiment, phase modulator 26 may be represented based 
on an optical power of light emitted by light source 22 and an operating phase bias 

10 of FOG 10. An operating phase bias can refer to a phase bias applied to 
counterpropagating light beams 34, 36 in optical waveguide 28 to displace the 
operating point of FOG 10. In one embodiment, the phase modulator 26 may be 
represented based on a product of the optical power and a sinusoidal function of 
the operating phase bias. For example, the phase modulator may be based on the 

15 product 

(Eq.2) K^sin ( <t> b ) , 

where I G is the optical power of light source 22 and <)) b is the operating phase bias 
of FOG 10. 

Feedforward component 130 may also include a representation 134 of a 
20 signal digitizer 30. Generally, a signal digitizer 30 may include a light detector, an 
analog-to-digital converter (ADC) , filter(s), and other processing component(s). 
A variety of signal digitizers may be represented based on schemes described 
herein. 

In one embodiment, the signal digitizer 30 may be represented as including 
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a photodetector and an associated preamplifier 135, a filter 136, an ADC 137, and 
a sampler 138. The photodetector and associated preamplifier 135 may be 
represented based on a photodetector scale factor R d , a preamplifier impedance Rf 
and a preamplifier gain G e . The photodetector scale factor R d may represent a 
5 scale factor between an input optical power and an output analog signal, e.g. 
current or voltage. In one embodiment, the photodetector and associated 
preamplifier 135 may be represented based on the product of the photodetector 
scale factor R d , the preamplifier impedance R^, and the preamplifier gain G e . The 
ADC 137 may be represented as a gain based on a number of bits b in the ADC 

10 137. In one embodiment, the ADC 137 may be represented as a gain based on the 
power 2 5 " 1 . In one embodiment, the filter 136 may be represented as a gain G f in 
voltage after the photodetector and associated preamplifier 135 and before the 
ADC 137. The sampler 138 may be represented as a sampler for analog-to-digital 
conversion. Accordingly, in one embodiment, the signal digitizer 30 may be 

15 represented based on the product 

(Eq.3) R - R l G.G f 2 w - 

Feedback component 140 may include representations of at least one FOG 
optical component and at least one FOG electrical component. Feedback 
component 140 may include a representation 142 of a demodulator 32. Generally, 

20 a demodulator 32 may include a sampler, a truncator, a digital-to-analog converter 
(DAC), and other processing component(s). A variety of demodulators may be 
represented based on schemes described herein. 

In one embodiment, the demodulator 32 may be represented as including a 
sampler 143, a truncator 144, and a DAC 145. The sampler 143 may be 
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represented as a sampler for digital-to-analog conversion. The truncator 144 may 
be represented as a digital truncation gain G D that occurs after the sampler 143 and 
before the DAC 145. In one embodiment, the digital truncation gain G D may be 
based on the number of bits d' in the sampler 143 and the number of bits d in the 
5 DAC 145. For example, the digital truncation gain G D may be based on the power 
2 d d . The DAC 145 may be represented as a gain based on the number of bits d in 
the DAC 145. In one embodiment, the DAC 145 may be represented as a gain 
based on the power 2 2_d . Accordingly, in one embodiment, the demodulator 30 
may be represented based on the product 

10 (Eq. 4) 2 d " d '2 2 " d =2 2 " d *. 

Feedback component 140 may include a representation 146 of a phase 
modulator 26. In one embodiment, phase modulator 26 may be represented based 
on a phase modulator scale factor Kp m . The phase modulator scale factor Kp m may 
represent a scale factor between an input analog signal, e.g. current or voltage, and 

15 an output angular measure. 

Feedback component 140 may also include a representation 148 of an 
optical waveguide 28. In one embodiment, the optical waveguide 28 may be 
represented as a time delay. The optical waveguide 28 may be represented as a 
transit time x for light to 

20 propagate through optical waveguide 28. For example, as previously described, an 
optical waveguide 28 may include a coil of optical fiber. In such an embodiment, 
the optical waveguide 28 may be represented based on a transit time 
(Eq. 5) x = nL / c , 

where L is the length of the coil and n is the index of refraction of the optical fiber. 
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Fig. 2 is a block diagram of an embodiment of an exemplary feedforward 
component for a closed-loop transfer function 100 according to Fig. 1. As shown 
in Fig. 2, feed forward component 200 may include disturbances at summing 
points 202, 204, 206, and 208 based on an optical power noise ^ 205, a shot noise 
5 i s 210, a preamplifier current noise i n 220, a preamplifier thermal noise i R 230, a 
preamplifier voltage noise i v 240, and an ADC quantization noise n^ 250. As 
described herein, a transfer function 100 having a feedforward component 200 
may be used to compute a coefficient of random walk (CRW) and a power spectral 
density (PSD) of noise of FOG 10. 
10 A PSD of shot noise i s 210 may be represented based on a photodetector 

current i D . In one embodiment, a PSD of shot noise i s 210 may be represented 
based on the product 

(Eq. 6) 2qi D =2qI 0 R D (1 +cos ((j) b ) ) , 

where 1°' R d , and (j) b have been previously defined, and q is the charge of the 
15 electron. 

A PSD of preamplifier thermal noise i R 230 may be represented based on a 
temperature T K of the FOG 10 and a preamplifier impedance Rf. In one 
embodiment, a PSD of thermal noise i R 230 may be represented based on the 
product 

20 (Eq. 7) 4 kT^R,, 

where k is Boltzmann's constant. 

A PSD of preamplifier voltage noise i v 240 may be represented based on a 
preamplifier voltage e n , a preamplifier noise gain G n , and a preamplifier impedance 
Rf.. In one embodiment, a PSD of preamplifier voltage noise i v 240 may be 
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represented based on the product 

(Eq. 8) (fifiJRf) 2 . 

A PSD of ADC quantization noise n ADC 250 may be represented based on an 
ADC sample period t, a preamplifier impedance R, a filter gain G f , and a number 
5 of bits b in ADC 137. In one embodiment, a PSD of ADC quantization noise n^ 
250 may be represented based on the product 
(Eq. 9) 2t/[12(R f G f 2 b " 1 ) 2 ] . 

PSDs of optical power noise I n 205 and preamplifier current noise i n 220 
may be represented based on schemes familiar to those of ordinary skill in the art. 

10 Generally, transfer function 100 may be manipulated using well known 

control system transform theory to determine relationships between FOG 
performance parameters and physical parameter(s) associated with FOG 
component(s). Appendices 1-5 include features related to manipulation of transfer 
function 100. Relationships for an operating frequency, a bandwidth, a PSD of 

15 noise, and a CRW are provided immediately below. As shown, these relationships 
may depend on FOG physical parameter(s) including at least one of an optical 
power I 0 of light transmitted by a light source 22, an operating phase bias <J) b , a 
photodetector scale factor R d , a preamplifier impedance R f , a preamplifier gain G e , 
a filter gain G f , an ADC gain 2 b "\ a phase modulator scale factor Kp m , and a transit 

20 time x. 

Based on a transfer function 100 having a feedforward component 130, an 
operating frequency co o for a FOG 10 may be expressed as 
(Eq. 10) co o = I Q sin(<t) b ) R d R f G e G f -2 b l • 1/x *G D -2^ ♦K pm 

Based on a transfer function 100 having a feedforward component 130, a 
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90° bandwidth BW90 for a FOG 10 may be expressed as 

(Eq. 1 1) BW90= ( 1 80/ti • arg (H(e Hp \ / 0 ))+90f , 

where H(z, I D ) is defined by the equation 

(Eq-12) H(z, I 0 ) = a> 0 .f ^ N+ " M+l) -z- 2 , 

5 in which N and M are described in Appendices 1-5, as those of ordinary 

skill in the art will understand. 

Based on a transfer function 100 having a feedforward component 130, a 3 
dB bandwidth BW3 for a FOG 10 may be expressed as 

(Eq. 13) BW3 = ^{H(e'™J 0 ))--j=j. 

10 Based on a transfer function 100 having a feedforward component 230, a 

PSD of noise for a FOG 10 may be expressed as 
(Eq. 14) 


PSD= 1 


{k s k,kJ 


2ql D + 


•2 Gl-el 1 2-t 

+ l n + 7 +/ V 

^ (ft r a f -2^y 12 


Based on a transfer function 100 having a feedforward component 230, a 
15 CRW for a FOG 10 may be expressed as 


(Eq. 15) Ci?W " ^ ' PSD 


n V 2 


2 

Based on the relationships provided in Eqs. 10-15, performance parameters 
for a FOG 10 may be computed. Generally, a performance parameter may be 
computed by substituting values of physical parameter(s) in the corresponding 
20 relationship for the performance parameter. For example, an operating frequency 
of a pre-existing FOG may be computed by substituting the values of the physical 
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parameters of the FOG in the relationship for the operating frequency provided 
herein. As previously indicated, physical parameters can include, for example, at 
least one of an optical power I c of light transmitted by a light source 22, an 
operating phase bias <{> b , a photodetector scale factor R d , a preamplifier impedance 
5 Rp a preamplifer gain G e , a filter gain G f , an ADC gain 2 b ~ 1 , a digital truncation 
gain G D , a DAC gain 2 2 ~ d , a phase modulator scale factor Kp m? and a transit time x. 

The relationships provided in Eqs. 10-15 may be used to design a FOG • 
having desired performance parameter value(s). In one embodiment, performance 
parameter value(s) may be provided. Based on the relationship(s) corresponding 

10 to the performance parameter(s), value(s) associated with physical parameter(s) 

may be determined for which the computed performance parameter(s) will have or 
approach the performance parameter value(s). Initial value(s) associated with 
physical parameter(s) may also be provided. The performance parameter(s) may 
be computed based on the corresponding relationship(s) and the initial value(s). If 

15 a difference can be determined between the computed performance parameter(s) 
and the performance parameter value(s), then value(s) associated with physical 
parameter(s) may be iteratively adjusted, and the performance parameter(s) may 
be iteratively computed based on the iteratively adjusted value(s). For example, a 
desired value of an operating frequency may be provided, and values of physical 

20 parameter(s) may be determined for which a FOG will have the operating 

frequency value. Also for example, desired values of an operating frequency and 
a PSD of noise may be provided, and value(s) of physical parameters may be 
determined for which the operating frequency and the PSD of noise approach the 
desired values. Generally, the relationships provided in Eqs. 10-15 may be used 
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with regression schemes familiar to those of ordinary skill in the art. 

The schemes described herein are not limited to a particular hardware or 
software configuration; they may find applicability in many computing or 
processing environments. The schemes can be implemented in hardware or 
5 software, or in a combination of hardware and software. The schemes can be 
implemented in one or 

more computer programs, in which a computer program can be understood to 
include one or more processor-executable instructions. The computer program(s) 
can execute on one or more programmable processors, and can be stored on one or 

10 more storage media readable by the processor, including volatile and nonvolatile 
memory and/or storage elements. The programmable processor(s) can access one 
or more input devices to obtain input data and one or more output devices to 
communicate output data. 

The computer program(s) can be implemented in high level procedural or 

15 object oriented programming language to communicate with a computer system. 
The computer program(s) can also be implemented in assembly or machine 
language. The language can be compiled or interpreted. 

The computer program(s) can be stored on a storage medium or a device 
(e.g., compact disk (CD), digital video disk (DVD), magnetic disk, internal hard 

20 drive, external hard drive, random access memory (RAM), redundant array of 
independent disks (RAID), or memory stick) that is readable by a general or 
special purpose programmable computer for configuring and operating the 
computer when the storage medium or device is read by the computer to perform 
the schemes described herein. 
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While the schemes described herein have been particularly shown and 
described with reference to certain exemplary embodiments, those of ordinary 
skill in the art will understand that various changes may be made in the form and 
details of the schemes described herein without departing from the spirit and scope 
5 of the present disclosure. 

For example, transfer function 100 may be modified based on schemes 
described herein to compute performance parameters of FOGs including 
components and/or arrangements of components similar to or different than those 
of FOG 10 shown in Fig. 3. 

10 Those of ordinary skill in the art will recognize or be able to ascertain many 

equivalents to the exemplary embodiments described herein by using no more than 
routine experimentation. Such equivalents are intended to be encompassed by the 
scope of the present disclosure. Accordingly, the present disclosure is not to be 
limited to the embodiments described herein and can include practices other than 

15 those described, and is to be interpreted as broadly as allowed under prevailing 
law. 
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Computation of 3 dB and 90° Bandwldths from Digital Model 
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Ph(a,,Io):=-^.arg(H(e i<DT ,Io)) 



BEST AVAILABLE COPY 
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0,50..2000 


FOG Transfer Function - Magnitude 



Frequency - Hz 

Nominal Detector Power 
Half Detector Power 
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A STUDY OF GAIN DISTRIBUTION AND RANDOM WALK 


The following four rules will be used to study gain distribution: 


1. BACKING OUT A NOISE TERM 
THROUGH A GAIN BLOCK 


-n(t)- 


-*.<t: 



J L EQUIVALENT 

(tyG— 


J > 


Signals. Noise Sources, and Gains 


WW) 


2. For a modulation/demodulation block containing the 
modulation function M(t), use: 


M(t/ 


or M(t) = M(t) 


-1 


3. If n(t) is a white noise random process, then for the 
purpose of statictisal computations: 

n(t) = M(t>n(t) 

since both terms have the same statistics. 

4. If Sm^©) is the PSD of a random process n(t), and 
M(t) is a square wave modulation function of 
frequency ©o, then, y(t) = M(t)-n(t) has PSD: 


n=±1, ±3, 


I 

^ADcW 


Feedback 


ADC 


H(t) is the input rate 

Ks * 2 7t<L ' P j s the Sagnac scale factor 

Kj b Io-sm(4b) is the phase gain at the operating bias point, <j>t, 

Io is the optical power (1/2 peak) 

In(t) is the optical power noise 

Rd is the detector scale factor in A/W 

Rf Is the feedback resistor in the transimpedance amplifier 

Ge is the net voltage gain from the detector to the A/D input 

G n is the noise gain of the transimpedance amplifier 

is(t) is the shot current 

iR(t) is the feedback resistor thermal noise 

in(t) is the amplifier current noise 

en(t) is the amplifier voltage noise 

nadc(t) is the A/D quantization noise 

ADC = is the gain, in Isb/V, for a A/D with b bits. 
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Performing this procedure for the Sagnac phase gives 


x ,* ■ M (»^C0 jgffl iR(0 i n (0 M(t).pn(t>O n 
4>s(0 + — " + „ ^ + + + — ■ + 


nadc(t) 


Kj KrR D Kj.Rd K,-R D Ki-RD-Rf ' K I .R D .R r Q E .2 b " 1 
A factor M(t) indicates those noise sources that may not be white. Backing out alt the way to the input rate gives: 


M(tMn(t) is(t) i R (t) 
O(t) + — + __ + 


in(t) M(t)cn(t)Q n 

+ __ — — + • 


nadcO) 


Ks-Ki Ks-Ki-Rd Ks-Kj.Rd K^-Rd Ks-KrR D -R f Ks* Kj- Rd * Rf • Ge- 2**~ 1 


This expression can be used to compute the net sampled PSD of all noise sources. The PSD of the shot current is: 
Pis(o) ■ 2.q-i D — where i D = io-(i + cos(<|> b )).RD is the detector current 


rt -19 


and q := 1.602- 10 is the electron charge 

i 

4 k T ^ 

The resistor thermai noise has PSD: PiR(a>) where k := 1.380658- 10~ 23 is BoltzrnarVs constant. 


Rf hz 


Tk := 298 is the Kelvin temperature is 


The A/D quantization noise has PSD: PacdU) c where x is the sample period. 

With the following parameter values: 

Rf := 30- 10 3 R D := .55 b := 8 1q := 6-10~ 6 Ge := 49-3.6 Ge = 176.4 


in -0-10 i: -pr e,, := 32-10 G n := 1 

yhz yfhz 


c:=310 n:=1.45 


Coij and wayelenqtl) friers); 
Transit time: % := — 


L := 600 D 


59.2 
1000 


i 0 A C ,«~9 „ 2-7T-L D 
X := 845-10 Ks := 


X-c 


t = 2.9x 10 6 t 1 = 3.448x 10 5 


Modulation period: 
The net PSD, with units 


T := 2-t T = 5.8 x 10 6 T 1 = 1.724 x 10 5 


can be written: 


\secj hz 

qJoU *cps|jffi)-Rp + I— I + in + z + " [^- ■^^A *- 



2 12 


and from the PSD. the random walk coefficient is: mifoik - «0<— • i£i 
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ID = Io(l + cos(4»b))-RD is ■ V 2 q»D 

Ki « Io sin(4>b) 

[(*S(t)M(t)KiRD + is(t) + iR(t) + i n (t))Rf + e n (t)-G n ]G B -2 b ~ 1 + nad c(t) 
M(t) 2 = 1 


M(Q.K I -R D -R P GB.2 b - 1 L s (t) ■ ^^13^1^1 , , M(t)- nadc (t) 

[ Kj-Rd Kj-Ro-Rf Kj. RD . Rf . GE . 2 l>-l 


But M(t)-is(t) - i S (t) M(t) iR(t) ■ iR(t) M(t) in(t) = in(t) 

M(t) en(t) ■ e„(t) M(t). na d c (t) = nad c (t) 


M(t)Ki-R D Rf-GE-2 


b-1 


. m , iS(t) + iR(t) + in(t) en(t) G n nadc(t) 


■) 


^SW + + + 1 + - 


Ki-Rd Kj.Rd I Kj Rd Ki 


( in(9 | enffl^Gn^ l 
^KjRd * Kj-Ro-Rf J 


nadc(t) 


KiR D .R r G E -2 


b-1 


+S(t) ° Ks-0(t) Ks * 


2*7tL»D 
Xc 


ks. oo + -S» + *» + («L + *£2l 1 + 

Ki-R D Ki-R D I^Ki-Rd Ki-R D RfJ 


nadc(t) 


R D Kj-Ro-Rfj Kl . RD . RrGE . 2 b-l 


Ks- 


Q(t) + _J^ + _ML_ + 


[ks-IQ.Rd + Ks-Ki-Rp-Rf J " 


nadcO) 


Ks-Ki.R D -RfG E -2 


b-1 


Pis(oJ * 2-q iD — 
hz 


PiR(G>) = — — 

Rf hz 


iD 55 Io (l +cos(<}> b )) RD 
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2- g-ip 
PO = + • 


& Rn) 2 I Rf J (Ks-KiB 


« 2 2 
G n e,, 


(Ks-IQRD) (Ks-KtRd) ^ Kf ^ (Ks-KiR D )^ (Ks-Kj-Ro-RfJ^ (j^r^ q^*-!) 


2t 
2 12 


losing) iD = Io .( 1 + cos (^)) 


(Ks-Ki-Rd) 


2 2 
•en 


1 2-T 


2qiD+ — +in + ; , 


*0' 


(Ks-Iosin(* b )R D ) 


2qIo(l + cos(* b ))-R D + 


4 ~k -T£\ . 2 G n 2 en 2 1 
+ in + r — + ■ 


2-T 

2" 12 


4> b :=0,.l..n 


* J 2.(Ks-Io sfa(4» b )-R D ) 2 


Rtheim(4>b) : = 60- 


180 


* J 2-(Ks-Io-sin(* b )-R D ) 


\ Rf J r. • ■> 


ainp_i(*b) := 60 — • 

71 ' 2-(Ks-K)-sin(*b)-RD) 


amp_v($bj : = 60 


n 2 2 
G n en 


75 J 2-(Ks-Io sin(* b ).R D ) 2 Rf 2 
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0.01 


RANDOM WALK vs PHASE BIAS 



• 


! / 








I • 
1 • 










\ \ *• 


* / I 

! •IT 





•* >v / 
•* if 1 

\ \ \ 


/f f 

' // 









\ 

f 
i 





1 10 


-3 


l-icr 


Total Random Walk Coefficient 

Photon Shot Noise 
A/D Quantization Noise 
Resistor Thermal Noise 
Amplifier Current Noise 
Amplifier Voltage Noise 


Phase Bias - (units of PI) 


40.10- 6 ^-i-^.3600.10- 3 = 0.01f^i.^L 
2 Ks it hr uA 
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Point Sensitivities for Non random Inputs 


■ n 


10 9 180 


Ks-(lo-sin($ b )) n 


3600 = 39.048 


deg I 
hr 'tjW 


— 1N @ Detector input 


10 


-9 


180 


Ks-(losin(((> b )).R D 


3600 = 70.996 


deg 1 
hr t|A 


1 N @ Detector output 


10 


180 


Ks- (io* sin(^b)) ' RD'Rf * 


3600 = 2367 


deg 1 
hr 'pV 


1 N @ Trans-impedence Amp. output 


10 


-6 


180 


Ks (lo* sin(*b)) R D* Rf • Ge 


3600 = 0.013 


deg 1 
hr "pV 


1N@A/D input 


V 0 := 10 6 n := 2 PPM := 100- 10 6 (tuning error) Kp m := — 

4 


2 Kpm 
x — - — 


PPM ._2 v 0 f — .360o) = 0.000484 

n 2 -l V* J fr'HV 


2N @ IOC input 


Herv6 Lefevre's A/D bits criterion: 


<D! := 2.71.400-10 3 <d 2 := 2-H-800-10 3 B L := ffilC>2 r = 4 i 8 . 9 KHz 

4-^0)!+ ©2) 1000 


log 2 (x) := 


log(x) 
log(2) 


Criterion is rms-noise - Isb /i?BL-Rf Ge = -j^y 


(based on shot noise only, 
since it should dominate) 


or b 


:= 1 + log2 1 I 

L Rf GE'V BL-2-q-l0-(l + cos(4> b )).R D J 


b = 9.1 bits 


hjamp ' = 19 bade - = 8 


« 1 (2^-1) 

tWp-1 Ks t 


180 dee , 

= 34155.7 — (max angular acceleration) 

sec 2 


b 0ll t - 6 


n 1 180 

wT'ks' T '~' 36 °° = °' 066703 810860 < LSB va,ue for b out output bits from ramp) 
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Rf-GE = 5.292 x 10 6 Iq = 6x 10 6 4> b = 0.5tc Ks = 0.88 


One Isb, at A/D b := 8 


180 


K S Iosin(* b ).RDRrGE-2 b " 1 n 


3600 = 104.811 


m 


bit 


t = 2.9x10 6 DAC:=™ *dig ; = ~ 

2 12 ^ 12 


128 


at 1st Integrator 


KpnrDA OKdig 180 (^\ 

— ~ ^3600«— = 2.808 I hr J 

*S n \ ./ 

bit 


Ks Io-sin((}> b ).R D .R r GE^ (b ^ n - =9.541x10" 3 tt 
v ' 180-3600 


10 


IO-sin($ b ).RD-RfGE-2 


= 447.359 
VI bit 


>-l 


IOSin(<|) b ) >R D .R r GE-2 _ 3 bit 
— = 2.235 x 10 * 


10 


jiRad 


10" 


RD'Rf'GE-2 
i 


b-1 


= 2.684 x 10 


r 3 jiWatt 
bit 


Rp-Rf -GE^ 1 bit 
= 372.557 


10 


jxWatt 


shot(lo) := RrGB-^BL-2-q-I6-(l + cos(4> b ])R D 

N:=100 i:=l..N I0i lo~ Vnnsjihoti := shot(lOi) Vbit := 


A/D Dither with Shot Noise 



1 10 


110 7 i-uf 
DETECTOR POWER - MICROWATTS 


110 
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FOG ANGLE RANDOM WALK 


P(a>) = P 0 
(PSD) 


25o 


White Nois 

F(s) 



A 



X 


A 



► 


s V P 0* B L <J S - ^Po-BL + CTquant 2 
Bi = noise equivalent bandwidth of filter F(s) 
2 Isb 

Equant 88 — (A/D quantization noise) 


CT I - T^N.(p 0 -B L + a quant 2 ) 


N = = number of samples in sum. 

x 

For N«x s T aV e = 1-hr. this gives the ARW coefficient 
(This formula assumes that the individual samples are 
uncorrected, which will be the case if the filter time 
constant is much less than the sample period.) 


Some Noise Equivalent Ban dwidths 
F(s) Bl (Hz) 

©0 


S + fflQ 

©1 ©2 


s + ©j s + ©2 
2 

©0 

~2 — - 

s +2-^-©o*s + ©o 


1-e 


-sT 


s-T 


©0 


©j-©2 


4-(©j +© 2 ) 

71 


i 

2-T 


(average over period T) 


Definition 


B L - J"(|P(2.w-i-f)|) 2 df 
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Setoaram ter values 
The PSD of the shot current is: 

A 2 


Pis(«>) B 2-q-iD — wher i D = ^-(l + cos(<|> b )).R D is the detector current 


and q:= 1.602 10 19 is the electron charge 


The resistor thermal noise has PSD: PiR (a>) = ^ where k := 1.380658- 10~ 23 is Boltzman's constant 

T K := 298 is the Kelvin temperature is Tk := 0 (set T K = 0 when thermal noise is included in 

With the following parameter values: ^ = - amplifier output noise spec) 

°' 2 

3 

Rf:=3010 R D :=.55 b := 8 Io:=610~ 6 G E := 49 Gp := 3.6 

. _ n , n - 12 A _ 9 v „ o (set % = 0 when cunrent noise is 

* ° 7^ Cn:=3210 7^ G n ;=1 c:=3-10 n := 1.45 included in amplifier output noise 

spec) 

Coil and wavelenoth (meters!: L := 600 D := — X 845-10"" 9 Kr - 2 ' 7t ' L ' p 

1000 ' ^ " x . c 

Transit time: <r := — x = 2.9xl0~ 6 -T 1 = 3.448 x 10 5 


c 


Modulation period: T := 2«x T=5.8x10~ 6 T~ 1 - 1.724 x 10 5 m := - f •= - 

T s ' 


T 


© M := 2-TC-Qvf © s := 2-7t f s 

i 2 2 

Th excess noise h?s PSD: Pcx ( a) ) = ^ where Af is the optical spectrum frequency width 

In terms of the full width at half maximum, Aff^^ Af is computed as: 

AX^:, 18-10-* Af^.^E M:= rZZ Mfwhm rjZ m 15Q5 

X V 2-ln(2) ^ 21n(2) 

Low pass filter 

• , := 2.,c.4.42.10 6 co 2 := 2-n-3.62.10 6 P(s) := B L := -fl^ 2 b L = 3.126* !0* Hz 

s + ©i S + CD2 4 ^©1+02) 

A/D bits: b:=8 gain ADC := — — (+/-2v range) 

2 v 


arcsec 


2«7i 1 180 

_ T 3600 = 0.066703 — 

2 6 n Isb 
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Compute ARW 


volt 


1. The net PSD before the filter F(s), with units of can be written: 

hz 


PSDo(lO>*b):= 2-q 
P 0 := PSDo(lo,*b) 


I 0 -(l+cos((( >b )).RD + L — ' + [~~ ^~ J +ln + _ 


RfW 


P0 = 6.81 x 10 


- 12 volt 


hz 


2. The rms filtered noise after F(s) is: 07 := Gf-^Pq-Bl o-f = 0.0166 vrms 


3. The rms sampled noise is: 


as = 2.146 Isb rms 


4. Accumulate samples for one hour and multiply by dt to convert to angle: 


180 1 / 3600 2 


7C (Ks-Io-sin(<>b)*RDRrGE*GF-ADC) 


ai = 0.006383 


dcg 


b^c := 8 AVh := 1 


^badc-1 
AVh 


•Gp . 


2q 


IO-(l +cos(* b ))R D + L — * 


2-22 
- + — 


Rf 


Rf 2 G E 2 B L = 2.126 


8.933 0.482 6.309 


0.871 


180 


n (Ks-Io-sm(<|> b ).RD) 


.>/3600-t%/Bl- 


2-q 


l0-(l + cos(ct> b ))RD + L - i •^ 2U ' 


l2 2 2 
G n -Cn 


Rf 


= 0.006325 


8.933 0.482-6.309-0.871. (0.000267) = 0.006317 
RD 1.226 Iq = 6x10" 6 
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Define functions for graphing 

CD l'©2 


. 180 [36Q0 /(ADC -Gp) 2 - 2-q Io (l +cos(^ b ))-R D R f 2 -GB 2 B L (a) i ,co 2 ) 

shot(lo,4>b.«M »a>2) == T- / ~i *7 — r~ — \ 

V 7 « V * (KsIo-6in(<t>b)RDRfGEGFADC) 


excess 


/ r , \ 180 13600 


1 


(ADCGp) 2 


[lb'(l+coe(» b ))-RD 


Af 


2 2 
•Rf 'Ge 


Bl(cdi ,co 2 ) 


(Ks-Io- sin^^RD-Rf-GEGF- ADC) 


/ \ 180 / 3600 


(ADCGp) 2 


Rf -Ge 


Rf 


Bl(cdi,co 2 ) 


(Ks-Io-sin(^b)-RDRf'GEGpADC) 


fj , x 180 / 3600 > I Rf J 

-P(l6.4b.»l.«2) := (Ks.Io-sin(H) RDRrG E GFAic) 


<t> b :=0,.l..it 


[3600 


(KsIosm(* b )R D RfGEGFADC) 


or lsb( I 0.<l>b> a >l «">2) == ADC-Gp-^PSDo^.^bj-BL^l ,«>2) lsb ms at in P ut to ^ 

180 nf J(ADC Gp) 2 -PSDo(lo ,» b ) B L (<D 1> a, 2 ) + -^ 
o»(lo.*b.-l.«2.T) := V- T VT (Ks-Io sin^^ RD Rf Gb-Gf-ADC) 
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— =4.42xl0 6 . — = 3.62xl0 6 In = 6x 10~. 6 
2-ji 2n 


RANDOM WALK vs PHASE BIAS 


\ 


; / 




y i 






Le, — . J 

H ! / 



i \c 



G B c ^\ 1 













• Total Random Walk Coefficient 
-* Photon Shot Noise 

* A/D Quantization Noise 
— Amplifier Noise 

■+ Excess Noise 


Phase Bias - (units of PI) 
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RANDOM WALK vs BREAK FREQUENCY 


CD i - ©2 








Iq = 6 x ] 




















; 

! 













0.007 
0.006 
ta 0.005 

i 

"° 0.004 
0.003 
0.002 
0.001 


100 


^ 10 


1 


Frequency - MHz 
A/D INPUT NOISE vs BREAK FREQUENCY 



2 2.5 3 

Frequency - MHz 

TIME CONSTANTS TO SAMPLE POINT 


3.5 


1 


1.5 


2 2.5 3 

Frequency - MHz 

Note: boxes are at multiples of the modulation frequency. 


3.5 


4.5 



4.5 
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Discussion of aliasing 

np := 100 i := 0..np fi := 10 np MFa(f) := |F(2-7i-i-f)| MFsai := MFa(fi) 


100 

£ (|F(2. 7 ci.n.f s )|) 2 =4.26 
n = -100 


F(s) Magnitude Transfer Function 
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0 0.340.691.031.381.722.072.412.76 3.1 3.453.794.144.484.835.175.525.866.216.55 6.9 7.24X597.938.288.628.979.319.66 1010.34 


Frequency - MHz 
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Optimization of SNR 


Step function for F(s)with ©i ■ ©2 s3 ©0 s(t,© 0 ) := l-(l -hct> 0 «t) e 


©0*t 


Signai-to-noise ratio: snr(ft) 


:= 10 4 ,2-10 4 .. 10 6 


Iq = 6x 10' 


-6 


SNR vs BREAK FREQUENCY 


\ / 

0" 

5 

2-(fM 

io- 6 )! 


J 

— 1 






















1.5 


0.5 


0.1 0.2 0.3 


0.4 0.5 0.6 

Frequency - MHz 


0.7 0.8 


0.9 
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ARW vs I n Plot (semilog plot) 

Set phase and break frequencies for plot <[> b : = ^ 

Set maximum Iq for plot: Iomax := 1 01 0~ 6 watts 


ARW vs Modulated Power 



10 (modulated power) in micro-watts 

" Total Random Walk Coefficient 
M * M Photon Shot Noise 
oaa A/D Quantization Noise 
° Amplifier Noise 
Excess Noise 


NGC-162/ 000388-280 

CD! := 2«7t-10 6 G> 2 •= 2-ti-IO 6 

np := 50 i := 1 ..np 10] := — -Iomax 

np 
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ARW vs to Plot floo-loa plot) 

Set phase and break frequencies for plot <t>b : = ~ ® i *= 10 6 a> 2 •'= 2-ic- 10 6 

Set mimimum power level and number of power decades for plot: Iomin := .1 • 10" 6 n_decades := 3 

— n_decades 
np := 10 i := 0..np I0i := Iomin* 10 np 


ARW vs Modulated Power 



10 (modulated power) in micro-watts 

Total Random Walk Coefficient 

K >( * Photon Shot Noise 
B-B " B A/D Quantization Noise 

Amplifier Noise 
' ' 1 Excess Noise 
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©1 := 2-n-MO 6 ©2 2-k-MO 6 Ii := 10~ 6 <|> b := 0,.l ..n 



Phase Bias - (units of PI) 


Solid curves have (counting from top curve) Io = 1. 2, 3, 4, 5, 6 t 7, 8, 9, 10 pwatt 

Dotted curves have (counting from top ) Iq = 12, 14, 16, 18, 20, 30 uwatt 


BEST AVAE LADLE COPY 
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-6 


<D 1 ?=2-7t ^ a>2:=a>i Io:=610 (f) b := 0,-1 ..n 


RANDOM WALK vs PHASE BIAS 



Phase Bias - (units of PI) 


Solid curves have (counting from bottom curve) © x = m 2 = 2, 4, 6, 8, 10, 12, 14, 16, 18, 20 x © M 
where co^is the modulation frequency 


BEST AVAILABLE COPY 
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ANTI-ALIASING FILTER BANDWIDTH TRIM TABLE 

Set modulation bias: <t>h := — Set desired rms noise at AID: 


rmsjioise := 0.8 Isb 


Set lower limit for filter comer freque ncy as multiple of modulation frequency: i ower := 2 

(For larger values of Io the computed value of the filter comer frequency may be too low, according to some criterion other 
than rms noise, like settling time etc. Setting this parameter defines the lower limit) 


Compute table 
ii:= 1..24 


I f := l ° 6 l3 | f - set ij :== root(ai sb (ir0.5-10" 6 ,<t) b ,2.7;.f,2.7t.f) -rms_noise,f) 


„-6 


OUT = 





1.018 

l0 ^corner 


4 

0.784 

4.5 

0.719 

5 

0.662 

5.5 

0.61 

6 

0.564 

6.5 

0.522 

7 

0.485 

7.5 

0.449 

8 

0.418 

8.5 

0.39 

9 

0.364 

9.5 

0.345 

10 

0.345 

10.5 

0.345 

11 

0.345 

11.5 

0.345 

, 12 

0.345 j 


Io is the modulated detector power 

^corner is the comer frequency of a two stage filter with the same corner 
frequency in each stage 


<^lsb[l0»*b»2-(2.7c.^ 1 ) > 2.(2.n.4i)] = 0.626 
^lsb^Y»*b>2(2n.^),2.(2.7i^)J = 0.487 
<»l8b[l6.*b. 3 -(2-it-ftf) f 3.(2-w-fi4)] » 0.766 


2jM 
10 6 


= 0.345 
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HARMONICS OF AN INPERFECT MODULATION SIGNAL 


MODULATION SIGNAL 


r = rising slope 
f = falling slope 



y = r*(t-T) 


Harmonics of the modulation signal shown 
in the figure will be derived. The signal 
differs from a perfect square wave 
modulation in three ways: 

1. Finite, unequal rise and fall times 

2. Duty_cycle * 50% 

3. Unequal high and low magnitudes 


y = -L 


Use the following transform-inverse pain 
rT 


1 f °° 
h(n)= T J ^O ^O ^Odt F(t)= h(n) exp(-in.<D 0 .t) 


Examples: 

rT 


™ exp(i n.©o t).sin(© 0 t) dt = ^ for (n = = 0 otherwise 


> 0 := 2 n. 


1 r T 

— • exp(vn-(D 0 -t)-cos(© 0 .t) dt ■ - for (n = +/-1), = 0 otherwise 


— • cxp(i n ©o*t) exp(i m-a>o t) dt = 6(n,-m) 
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The transform of the modulation function shown in the figure is: 
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H 

HARn(T,tc,r,f ,H,L,n) = ~ exp(ina>o-t)(rt)dt ... 


0 

H 


c 


exp(in«a>o-t)(H) dt ... 


r 

L 


T „ ""l"-' 
•'t- — 


+ -• exp(ina> 0 t)[-f(t-tc)]dt... 
'«c- 


f 

,T-i 
r 


+ Y cxp(i-n.© 0 -t)-(-L)dt . 


L 

l 


**7 


+ ^- exp(i n © 0 t) [r(t - T)] dt 

JtJ: 


Carrying out the integral gives the general expression: 


HARn(T,tc,r,f,H,L,n) := »- 

4-n -n 


For n=0, this is: HAfa(T,1c.r,f ,H,L,0) « — ^ ]3±±H + (L + H)--L 

v 7 2 fr-T ' T 


rT-£exp£-2i-n«7P 
+ f-T- exp 2 i n-7i 


(-Tt + L) 
/ (Tr) 
CT-f) 


- exp 


- exp 


H Tl 
2-i-n-n 

L <? r >J-J s-,-, 

K- f+H ) 11 


-2i-nn- 


(Tf) 
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For r, f -> oo: har* 


r. nn o 

i_ir(T,tc,H,L,n) := _L.(L + H ). e V T Asin| n « - ] 
v ' nn \ 1) 


For n = 0. this is: 


(H + L) 


tc-L 


which gives the Fourier representation: 


M(t) 


• ( M 

* oo sm n-n — I - 

= (L + H)-|-L + ^.(L + H). J] -A_^. co J n . ao r t _i. 1 jl 

n = 1 


Note, that for ■ -T (50% duty cycle), this reduces to: 

mi 
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Onstantaneous rise and fall) 


M(t) = — -+--(L + H> V suifn-con-t) 
2 n ^-f n v y 

n= 1 

(n odd) 


(for the case where L * H is the only 
defect there are no even harmonics) 
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EVEN HARMONICS FROM DUTY CVuLE VARIATION 

If L = H = A in addition to istantaneous rise and fall, the representation becomes: 
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(Duty_cycle * 50% is the only defect 
here) 


or 


Ms(t) ■ 2- A 


c 4-A ^ sin(n-7c D) 

•o + y 

n ^-f n 
n= 1 


cos n-a>o 


{-HI 


tc = T-D D = ~ + 8 

2 


1 


6 ■ D - - (D is the duty cycle, D = 0.5 for 50%) 


where 

Picking out the even harmonics: 

MScven(t) = 2A-8 + ^ £ sip ( 2 ^' 5 ) -cosfr 1kd 0 t - 2 k^-5) 
71 k=l 

The 2nd harmonic term is: 

M$2(0 83 2-— •sin(2-7f5)-cos(2a> 0 t-2Tc*6) 

The even harmonic magnitude divided by fundamental magnitude is: 


4-A sin(n n-S) 


Relative_Magnitude(ii) ! 


sin(n-K-S) 


RM(6,k) ==20-loJ gE&Ottj) 


2 k- 


k:= 1..25 


(n even) 


66 k := 10 6 -l0 2S 


Relative Magnitude 


RM(.l,k) 
3 RM(.01,k) 
B RM(.001,k) -50 
S RM(.0001,k) 


-100 



i i 







1 

i 


Relative Magnitude of 2nd Harmonic 


0 20 40 60 

2 k 



M0 U 1 -10 


For 6 « 1: 


Harmonic Number (even only) 
M52W * 4.A-6-cos(2 <Do t-2 Tt-6) 


1 lO*" 4 1 10 3 
Duty Cycle (delta) 


0.01 


0.1 


Relative_Magnitude(n) = Tt-6 


01 = 2,4,6.... ) 


BEST AVA'LASLE COPY 
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WAVEFORMS 


Numerical study of even harmonics due to Duty_cycle * 50%: 

A := 1 D := .6 6 := D - — T := 1 a> 0 := 2-*-— N := 1000 

2 T 


i:=0..N tti := — T 
N 


100 . / x 

Ms(t) := 2 A-6 + > -cosj 

71 *-~t n 
n= 1 


[na>o-[t-i.D-T)j 


M6i := M5(tti) 


M5even(t) := 2-A-6 + S^±^ -cosfrk-CDo-t- Ik-**) 


k= 1 


MSevenj := M5even( tt i) 


MSodd(t) := 


4-A 


k=0 L v /. 


MSoddi := Maoddfa) 


D = 60% sinnal reconstructed from Fourier saries 
2 


M5j 0 
-1 

-2 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
tti 


MSoddj+MSeveiij 0 


-2 


Sinnal oonstmcted frnm even harmonics onlv 
1.5 


1 

MSevenj 0.5 
0 

-0.5 


! 


i 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
tti 

Sional constructed from odd harmonics onlv 
2 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


1 

18oddj 0 
-1 
-2 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
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EVEN HARMONICS DUE TO UNEQLmL RISE AND PALL TIMES 
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Setting L = H = A and = - T in the general formula, gives the harmonics for the case where th only defect in the 
modulation signal is unequal rise and fall times: 


HARn(T,r,f ,A,n) » — 
4 


r-exp 


(rT) J *1 (r-T) 


+ f*exp 


(f-T + 2 A) 
(f'T) 


(»v) 


- f-exp 


andforn = 0: HARn(T,r,f ,A,0) = 0 
This expression can also be written: 


HARn(T,r,f,A,n) = y - j rsinf MwA-- H-lf-sJlnn A— 1 

2-n L v T y V T / 

HASn(T,r,r,A,n) « T»i( 1 - rj. 


For r = f, this reduces to: 


Introducing = — tfeii = — the general formula becomes: 

This formula can be used to derive the following series for the even and odd harmonics separately: 

Expanding in the small parameters r" 1 and f " 1 (or equivaiently tnse , and tfeu ) gives the even harmonics: 


HARn(T ,t, f , A ,n_even) = in *' A 

6T 2 

which gives the representation: 


(2-A) 2 (2-A) 2 


i-n-rc-A 


(* 2 ♦ ^ 
tfell ~tee 


(small r 1 and f 1 
or small tee and tfau ) 


(f-T-2-A) 
(f'T) 


Meven(t) 


7C-A 


( * 2 ► 2 
tfall -tee 


T 2 ^ 


• ^ 2-k«sin(2'k-a>o-t) 


(small r 1 and f" 1 ) 
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WAVEFORMS 


T := 1 A := 1 ^ := .1 


tfeu := .2 <d 0 := 2 - 


- 7-|. j£?fe4»* + »- 3 ?]*=4»- 


N := 1000 i:=0..N tti := — T 

N 


:= Mevcn(tti) Moj := Modd(tti) 


Sinnal nonstnintftd from «vfin harmonics onlv 
2 


1 
o 
-1 
-2 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
tti 


2 
1 

Moi+Mej o 
-2 


I 


J 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
tti 


Sional oonstn ictp.ri frnm odrl harmonic*; onlv 
2 


1 

Moj o 
-1 
-2 


f 


j 


J 







0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
«5 
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HARMONICS OF A MODULATION 
SIGNAL 
WITH ASYMMETRIC DROOP 
MODULATION SIGNAL 



*» t 


y =-1 + b*(t - T/2) 


Use the following transform-inverse pair 

f 


Harmonics of the modulation signal shown 
in the figure will be derived. The signal 
differs from a perfect square wave 
modulation by linear droop shown in th 
figure, where the droop slopes are, in 
general, not equal: a * b . 


1 f T 00 
h(n) = -J exp(i. n .c>o't).F(t)dt F(t) = £ h(n)exp(-i.no> 0 .t) 


(on 2-n — 
T 


The transform of the modulation function shown in the figure is: 
T 

i f 1 

HARn(T,a,b,n) » -• exp(in©ot)o - *t)dt ... 
2 

Carrying out the integral gives the general expression: 

HA*n(T,a,b,n) = _2_.[[ 4 . n . _ J _ i. T .(a + b)-[(-l) n - l]] -n n.T.[a (-l) n -b]] 


4-n -n 


For even n this is: 


HARn(T,a,b,n_even) = — — T-(a-b) * — — -fa— -b--1 


t 
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which giv s the following Fourier series for the even harmonic part of the modulation function: 


n 

n=2. 4. ... 


If the nominal magnitude were +/- V 0i then the even harmonic magnitude relative to V 0 is: 

^(a.V 0 .l-b.V 0 ^l 

Relative Jtfagnitude(n) = -22Li ± V_ = J_ AV 

Vo nn V 0 

where AV = a-V 0 | - b-V 0 -| is the differential voltage droop. 

Using 201o^~J = -15.964 , the second harmonic can be expressed in dB as 

Relative Jvlagnitude_dB(2) ■ 20-logf — | - 15.964 dB 


Derivation of dosed fo rm approximation of ovrp bias error 

1 . The switch waveferm is: positive section negative section 

1-a-t 

2. Delaying with a tuning error gives: 


before delay: 


positive section 
1-a-t 


after delay: 


3. This is equivalent to: 


after delay: 


negative section 
-l+b(t-A) 


negative section 

-i + b(t-I.x) 

positive section 
l-a{t-I.T-A) 


positive section 
l-a(t-I.T-A) 


negative section 
-l+b(t-T-A) 
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4. Change sign of delayed signal: 

positive section 
1 - at 


before delay: 


after delay 
and sign change: 


positive section 
l-b-(t-A) 

5. Now add the two signals and multiply by tc/4: 
positive section 


f-f.(a + b)t + fb.A 
2 4 4 


negative section 
negative section 


negative section 


6. The positive section will be sampled at t ■ — T, and the negative section at t = --T 

8 " 8 


positive section 


negative section 


iS erenCe in det6Ct0r 9iveS 016 cosine of to* 18 * s '9 nals - and demodulation and integration takes the difference of the 


RECT - cosjj- - ^.(. + b) |.T + ib- aJ - cos^f + £.(a + b).|-T - J*a] 


8. Expanding in a and b: 


RECT* —.(a- b) 
4 


9. Applying the Sagnac scale factor and converting to deg/hn 


BIA S =I.r^. ( a-b)l.-^f-.I«2.3 6 00 
2 L 4 J 2-tc L-D 7i 


